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ABSTRACT 

We present a survey, using the Chandra X-ray observatory, of the central gravitating 
mass profiles in a sample of 10 galaxies, groups and clusters, spanning ~2 orders 
of magnitude in virial mass. We find the total mass distributions from '^0.2-lORe, 
where Re is the optical effective radius of the central galaxy, are remarkably similar to 
powerlaw density profiles. The negative logarithmic slope of the mass density profiles, 
a, systematically varies with Re, from a ~2, for systems with Re~4 kpc to a ~1.2 
for systems with 30 kpc. Departures from hydrostatic equilibrium are likely to be 
small and cannot easily explain this trend. We show that the conspiracy between the 
baryonic (Sersic) and dark matter (NFW/ Einasto) components required to maintain 
a powerlaw total mass distribution naturally predicts an anti-correlation between a 
and Re that is very close to what is observed. The systematic variation of a with Re 
implies a dark matter fraction within Re that varies systematically with the properties 
of the galaxy in such a manner as to reproduce, without fine tuning, the observed tilt 
of the fundamental plane. We speculate that establishing a nearly powerlaw total mass 
distribution is therefore a fundamental feature of galaxy formation and the primary 
factor which determines the tilt of the fundamental plane. 

Key words: Xrays: galaxies — galaxies: elliptical and lenticular, cD — galaxies: ISM — 
dark matter — galaxies: fundamental parameters 



1 INTRODUCTION 

X 

H ; 

. . 1 The global optical properties of giant elliptical galaxies are 
often parametrized by three key quantities, the effective 
(ha1f-liffVit) raHiiiH (R„) of the stellar light, the mean siir- 
face bri lghtness for luminositv) and the central line-of-sight 



must be tig htly correlated wit h the optical properties of the 



galaxy (e.g. Ciotti et al. 1996| ) 



Strong observational evidence for the existence of mas- 
sive dark matter halos around early-type galaxies has been 
provided by independent observational constrain ts from X- 
ray studies , lensing a,nd s t ellar dynamics (e.g. [Humphrey 



et al. 



2006, hereafter H06 



Mathews fc Brighcnti 



velocity dispersion (cro)- In this three-dimensional p arameter 
space, they occupy a narrow "fundamen tal" plane (Dressier 



et al. 



1987 



Djorgovski & Davis 1987), the orientation of 



references therein; [Kochanok _ 1995 ;_ Glrifhths 



2003 



ot al. 



and 



199f 



Gavazzi et al.||2007| ; [Gerhard et all^OOl^ [Thomas ct aljpOOTD 



which differs significantly from naive expectations from the 
virial theorem (assuming homology and a constant mass- 
to-light ratio, 7). The relative importance which deviations 
from homology, variations in the stellar population proper- 
ties and changing dark matter fractions play in producing 
this " tilt" lia,ve been ho tlv d ebated Ce.g. jRenzirii Cio tti 
1199^: Niorth ^ Ma,d;^B Ciotti et a.lJl199fl: IPrn..iel 



Iraham fc CoUes j 119971: 



Padmanabhan et al. 



et al 



p006| ; [Bohon et al.| p007| ; 



# 



2004; rrujillo et al 



Tortora et al 



Gerhard et al.[ ]2001| 
2004t PappeUari 



Although the total mass distribution within Ro is dominated 
by the stars, the dark matter fraction is non-negligible, and 
stellar dynamics and lensing studies (most of which are re- 
stricted to within ~Re) suggest that it establishes a conspir- 
acy with the luminous matter to produce total mass density 
(Pm) profiles close to pm oc (wher e R is the radius), i.e . 
the "singular isothe rmal sphere" (e.g. |< ochanck| 1995; Tret 

' ]2009( [Kronawit 



tor ct al 



Koopman; 

pooo^ 



2004; Koopmans et al 



200e 



200E, and 



references therein), with recent wor k tending to emphasize 



the importance of dark matter (e.g. Cappellari et al. 2006 
Bolton et al. 2007). To maintain the thinness of the funda- 



This "bulge-halo conspiracy" is similar to 
that which establishes flat rotation curves in disk galaxies. 

Total mass distributions approximately consistent with 
Pm OC have long been reported from hydrostatic X-ray 
analysis of n earby elliptical galaxies at sca l es much larger 
than Rc (e.g. iTrinchicri et al. 19861; [Thomasl 19861; pcrlcmit 



mental plane, however, the dark matter fraction within ~Rc 



SOS et al. 



199S; Nulsen fc Bohringer 199. 



Kim fc Fabbianc 
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1995; Rangarajan et aL||l99£| ; Matsushita et al. 1998) and, 
in their R osat imaRin g analysis of two galaxies, Buote & 
Canizai[e^ (1994, 1998) actually found an isothermal sphere 
potential to be preferred over some other mass distributions. 

With the imprr.^^crl fapgKilitioa nf Chnn^rn gnrl Yl\/rA/f miifti 



tighter honstraints on radial mass profiles have now been re- 



ported for a wide r array of galaxies (e 



200e| |2004 |2009at [Fukazawa et al.1^006 



po|04 P^Uivan et al 



3 'Sullivan & Pon 



2007: puote 



Humphrey et al 



2002 



larly resemble p,„ oc R ^ (as pointed out by Gavazzi et al 



which simi- 



2007 ), although not e xactly so over a l l rad ial scales ( |Ro- 



manowiiky ct al. 200t ). Fukazawa et al. (2006) fitted a model 
of the form pm oc R~" to the combined mass profiles of their 
galaxy sample outside 10 kpc, finding a = 1.67 ± 0.33. 

At the scale of massive galaxy clusters, hydrostatic X- 
ray analysis has revealed much l ess cuspy mass distributions 
(q ~1-1.4; e.g. [Buote fc Lewi^ ^OOj; [Voigt fc Fabian| |2006i 



Arabadiis et al.|pOO ^ avid et al. pooiT In particular, iLewid 



et al. ([^003 ) and Zappacosta et al. (2006) studied the very 
relaxed clusters A 2029 and A 2589, finding that the total 
mass profiles were in good agreement with the NFW shape 
expected for the dark matter only, once again suggesting 
some kind of conspiracy between the luminous and dark 
matter to produce an approximately powerlaw total mass 
distribution in the core of the system. Stellar dynamics and 
lensing studies of clusters have, similarly, found less cuspy 
total mass profiles than would be expected for an unmodi- 



fied NFW component pl us the stellar mass (e.g. Band et al 
2004 [Kelson et aL||2002|). 



In order to investigate this disparate behaviour at dif- 
ferent mass scales, in this paper we carry out a uniform 
X-ray analysis of a sample of relaxed galaxies, groups and 
clusters to investigate the shape of the innermost mass dis- 
tributions. X-ray analysis is ideally suited for this study, 
since it allows straightforward mass measurements in indi- 
vidual systems over a wide radial range, from the baryon 
dominated regime (R<Rc) to regions where the dark mat- 
ter dominates the gravitating mass (R3>Ro), thus providing 
adequate leverage to elucidate any conspiracy between these 
different components. For relaxed systems, hydrostatic equi- 
librium is believed to be an excellent approximation, with 
nonthermal effects contributing no more than --^20% of the 



total pressure (e.g. Churazov et al. 200^ ; Nagai et al. 2007 



after |H09|) . In our previous paper s 
^l[)07[ 



Piffaretti fc Valdarninj ^008[ ; [Humphrey et al.| |2009a|, here 



Gastaldello 



et al, pi|)07 ; Zappacosta et al. 2006) we have carried out de- 
tailed decompositions of the mass distribution of these sys- 
tems into baryonic and non-baryonic components. We have 
not, however, investigated in detail the relationship between 
these two components. In light of the apparent conspira- 
cies between them at both cluster and galaxy scales, in this 
present paper, we adopt the more pragmatic approach of 
examining whether the inner part of the mass profile can be 
modelled as a powerlaw, and investigating whether its slope 
varies systematically. 

Throughout this work, we assume a cosmology of 
Ho^70 km s"^ Mpc"\ ^^=0.3 and A=0.7. AU error-bars, 
unless stated otherwise, correspond to 1-a confidence re- 
gions. 



Table 1. Properties of the galaxies, groups and clusters in our 
sample. We list the K-band effective radius of the centra l galaxy 
(Rc) obtained from the 2MASS database (Jarrett |200C[ ), except 
for the clusters, for which w e report (t) the V-band Ro from 
Maluniuth fc Kirshnci ( 1985| ) and (*) the R-band Re from [Json[ 



et al\ { 1991 ). We quote luminosity distances for each object, based 



on the redshift in the NED database, except those marked "*", 
w hich were bas e d on the surface br ightness fluctuations study 



of ronry et al. (DOM, see [HOfij; H09). Deprojected density and 



temperature profiles were taken from the listed reference (Ref). 
Where the data are unpublished or were re-reduced in the present 
work, we list the Chandra observation ident ifiers (ObsID) and net 



exposu re tim e (Exp). Referenc e s: (1) H09, (2) Humphrey ct al, 

F 



priv. cor rim.) 
comm.); 



(2009b), (3) Humphrey et al.jMOSL 



Liu et al 



Gas t aldell o et al 
( [2009[ ) 



(4) F. Gastaldello (2009, 
(p007|), (5) W. Liu (2009, priv. 



Object 



Rc Distance 
(kpc) (Mpc) 



ObsID 



Exp 
(ks) 



Ref 



NGC1332 
NGC 720 

NGC 4649 
NGC 4261 
RXJ1159-I-5531 
MKW4 
AWM4 
ESO552-020 
ABELL2589 
ABELL 2029 



2.7 
3.1 

3.2 
3.4 
9.8 

10 

10 

16 

33t 

76* 



21.3* 
25.7* 

15.6* 

29.3* 

368 

87 

139 

138 

183 

350 



7062,7372 99 
8448,8449 



3234 
9423 
3206 
7190 
4977 



30 
71 
19 
52 
77 



(1) 
(2) 

(3) 
(1) 
(4) 



(5) 



2 DATA ANALYSIS 



2.1 The sample 



We chose a sample of 9 objects, spanning ~2 orders of mag- 
nitude in virial mass (Alvir) from the su rvey of morpholog l 



cally relaxed. X-ray luminous systems by Buote et al. (2007). 
We supplemented this sample with NGC 1332, which has the 
smallest Ro of the relaxed s ystems we have previously stud- 
ied for X-ray mass analysis (H06; H09|). We focused on Chan- 



dra data, since high spatial resolution is essential to study 
in detail the inner part of the X-ray halo. As we require 
coverage from < Rc to several times Rc , we only considered 
objects which had deep enough data to enable spectra to 
be obtained in at least two annuli within Rc. Unfortunately, 
most nearby galaxy clusters which are morphologically re- 
laxed on large scales also exhibit significant AGN-induced 
cavities in the core which complicate the analysis. There- 
fore, we only included two clusters which have been shown 



to be relaxed at both large and small scales ( Lewis et al 



2002; Zappacosta et al. 2006). The properties of the sample 
and the archival data we used are summarized in Table 



2.2 The Slope of the Mass Profile 

The Chandra data were reduced and analysed to provide 
de proje cted temperature and density profiles, as outlined 
For NGC 1332, NG C 4261 and NGC 4 6 49 the data 
and 



H09 



are described in 



HOE 



Humphrey et al. 



( 20081). We 

will d i scuss i n more detail the NGC 7 20 data in Humphrey 

For 



et al. 



(2009b) and the A 2029 data in 



Liu et al 



( [2009| ) 



RXJ1159-f 5531, we used the deprojected profiles described 
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Figure 1. Radial mass profiles for each object, arbitrarily scaled for clarity. The solid lines are the best-fitting profiles determined from 
our "forwa rd fi tting" analysis of the temperature and density profiles, while the data-points are determined from the more "traditional" 
method (§ We stress that the models are not fitted to these data-points but are derived independently. 



in Pastaldello et al.| ( POOT] , provided by F. Gastaldello 2009, 
priv comm.). 

Under the hydrostatic approximation, we transformed 
these density and temperature data into mass constraints 
by two complementary approaches. First, the "traditional" 
method involves parametrizing these profiles with arbi- 
trary models (for mor e details on these mode ls, see 



Humphrey et al. 2008; Gastaldello et al. 2007), which are 



H09 



then differentiated an d inserted into the equation of hydro- 
static equilibrium (e.g. Mathews 1978). An advantage of this 



method is that it makes no a priori assumption about the 
form of the mass distribution. By evaluating the resulting 
mass model at a number of radii (corresponding to each 
spectral extraction region) we obtained the mass "data- 
points" shown for each system in Fig 111 Error-bars were 



estimated via a Monte Carlo technique (iLewis et al. 2003) 



We here focus only on the central part of these data; based 
on experimentation, we considered the mass within lORo or 
200 kpc, whichever is smaller. Over this radial range, the 
profiles are all approximately powerlaw in form, but the ex- 
act slope varies from object to object. 

We find overall good agr e ement with previously pub- 
lished r aass profile s (|H0^; pI09|; [Gastaldello et alj|2007| ; ^ap 



pacosta et al 



2006 



Lewis et al. 



2003), although for ES0552- 



020 the normalization is 
by 



Gastaldello et al 



0.1 dex higher than that found 
using XMM data. 



( |2007[ ), using XMM data. Neverthe- 
less, this discrepancy is comparable with our estimated sys- 
tematic error for this object (§ 2.3), and will not affect our 
conclusions. 

Since the traditional method relies on the adoption of 



ad hoc temperature and density profiles, this can lead to sig- 
nific ant s ystematic errors in the recovered mass distribution 
(e.g. H09). Furthermore, the individual "data-points" are all 



correlated, which makes it difficult to interpret a fit made 
directly through them. Therefore, to be more quantitative, 
we fitted the mass dist ribution using t he "f orward fitting" 



method described in Humphrey et al. (2008). This involves 
solving the equation of hydrostatic equilibrium to compute 
temperature and density profile models, given parametrized 
mass and entropy profiles. Since the entropy profile must 
rise monotonically, we parametrized it as a constant plus a 
powerlaw with one or two breaks added, as needed. For the 
stellar plus dark mass distribution, we adopted a powerlaw, 
corresponding to pm oc R~°', i.e. 



M(< R) = M- 



R 



'"^ (rskpc 



(1) 



where M is the mass enclosed within radius R. M75 and a 
were parameters of the fit. An additional gas mass compo- 
nent was included self-consistently in the calculation, but 
is generally small in the fitted radial range. We fitted only 
the inner parts of the density and temperature profiles, as 
described above, freely varying logM75, a, the parameters 
describing the entropy profile and a te rm r elated to th e gas 



H09 



pressure at a suitable reference radius (HOE). Following 
parameter space was explored with a Bay esian Monte Carlo 
methodQ assuming fiat priors (see § 2.3 for the impact of 



^ Sne cificallv. the nested sampling algorithm of Feroz & Hobson 
(200S) 



4 P. J. Humphrey and D. A. Buote 



Table 2. Marginalized best-fitting values of the powerlaw slope 
(a) and normalization (Afrs) for each object in the sample. We list 
both the marginalized value and the l-tr statistical errors ("Best") 
and an estimate of the systematic uncertainty ("Sys.") due to 
various data-analysis choices (§ 2.3). We stress this should not be 
added in quadrature with the statistical errors. Also shown is the 
per degree of freedom (dof) for each fit to the temperature 
and density profiles. 

Object a logM^s 

Best Sys. Best 



Sys. 



NGC 1332 
NGC 720 
NGC 4649 
NGC 4261 
RXJ1159 
MKW4 
AWM4 
ES0552 
A 2589 
A 2029 



2.05 
1.91 
1.97 
2.01 
1.67 
1.66 
1.64 
1.39 
1.23 
1.21 



-1-0.05 
0.06 
-1-0.05 
0.07 
-t-0.03 
0.05 
-1-0.03 
0.05 
-1-0.11 
0.10 
-1-0.03 
0.04 
-1-0.13 
0.10 
-1-0.11 
0.18 
-1-0.17 
0.10 
-1-0.05 
0.06 



±0.04 
±0.02 

±0.03 

+0.09 
-0.04 

±0.06 
±0.03 
±0.08 

+0.09 
-0.07 
+0.27 
-0.10 
+0.22 
-0.11 



I2.l5i°;i? 



,2 34+0.08 
12 48+''-0* 

ir, q-, + 0.09 

^+0.10 
-0.05 
^+0.03 
*-0.02 

12.88lH^ 
12.8lin^ 

s+0.03 
-0.03 



12.82^ 
12.8 



12.i 

-,o ,7+0.01 
-^■'•-^'-0.01 



+0.17 
-0.07 
+0.03 
-0.03 
+0.23 
-0.09 
±0.02 

+0.07 

+0.06 
-0.02 
+0.11 
-0.06 

+0.10 
±0.12 
±0.11 



xVdof 



16/6 

5.7/8 

30/22 

6.5/2 

9.3/4 

26/18 

15/12 

9.6/10 

19/17 

23/8 




1 10 
Radius (kpc) 

Figure 3. Left: Comparison of the mass profiles of NGC 4261 
(circles), MKW4 (triangles) and A 2589 (stars), normalized to 
have the same mass at 100 kpc. Note that the shapes of the mass 
profiles are different even over overlapping radial ranges. Right: 
The same, but scaling the radius with respect to the virial radius 
(Rvir) and renormalizing each profile to have the same mass at 

O.OOlR^ir. 



using other priors), and the marginalized best-fitting mass 
profile parameters and l-cr confidence regions are given in 
Table y. Also shown is the per degree of freedom for the 
fits to the temperature and density profiles. In general, the 
fits are good. For the three objects for which the fits are 
not formally acceptable (null hypothesis probability <5%), 
the mean absolute fractional discrepancy between the tem- 
perature and density models and the data are less than 
~10%, indicating that the true mass distribution is, nev- 
ertheless, very close to a powerlaw in form. For one of these 
objects (NGC 4261), we did find larger discrepancies within 
~0.5 kpc, but we ignored the data within this region during 
our fit so as not to bias our results. (It is interesting to note 
that the full temperature and density profiles of this object 
can be fitted well by a mass deco mposition into luminous 
and dark matter components: 



H09 



For the two clusters, we 
have previously published powerlaw fits to the central mass 



distributi ons, based on other data ( Lewis et al. 



2003 



Buote 



& Lewj^ |2004{ [Zappacosta et sH. 2006 ) , and these studies 



agree well with our results. 

In Fig |l[ we overlay the best-fitting models onto the 
data-points obtained from the traditional method, finding 
overall good agreement. The mean absolute fractional dif- 
ference between the model and "data" varies from ~5% to 
~18%, and is smaller than 11% for 6 objects. While we don't 
expect a pure powerlaw to be an exact description of the 
mass profile, these residuals indicate that it is, nevertheless, 
an excellent approximation to better than ~10-20%. In fact, 
the discrepancies we found between the mass profiles deter- 
mined from both approaches are comparable to the typi- 
cal systematic uncerta intie s associated with the traditional 
analysis method (e.g. H09), so that the actual agreement 



between the mass distribution and a powerlaw shape may 
be even better. From Table ^, it is immediately clear that a 
systematically varies with mass scale. This is shown explic- 
itly in Fig ^, where we plot a versus Rc and versus M75 , in 
both cases revealing striking anti-correlations. 



2.3 Systematic Errors 

As in all studies of dark matter in early- type galaxies, re- 
gardless of the specific method, our work involved a number 
of arbitrary analysis choices. In this section, we describe how 
sensitive our results are to the various choices we made. For 
a more detailed discu ssion of these various systematic error 
assessments. 



H09 



We first examined how the choice of priors might be 
influencing our results by replacing each of the flat priors 
we adopted by priors which were flat in logarithmic space 
(we replaced the fiat prior on log M75 with one which is fiat 
on M75). These choices typically had a smaller effect than 
the statistical errors. Next we investigated the sensitivity of 
our spectral-fitting results to our treatment of the Chandra 
background by using the standard "background template" 
spectra, suitably renormalized to match the data at ener- 
gies > 10 keV, instead of the more robust, modelled back- 
ground adopted by default. To assess the importance of the 
radial range we fitted, we tried reducing it by ~20% for each 
object. We experimented with spectral-fitting over different 
energy ranges (0.4-7.0 keV, 0.5-2.0 keV and 0.7-7.0 keV, 
in addition to our default, 0.5-7.0 keV), varying the neutral 
galactic column density by 25%, and the distance by 30% 
(for the NGC objects, we instead varied the distance by the 
statistical e rror on their distance measurements in Tonry 
In Table we list the largest change in the 



et al 



2001) 



marginalized parameter values arising due to these choices. 
None of these systematic errors are large enough to affect 
our conclusions. 



3 DISCUSSION 

3.1 A luminous-dark matter conspiracy 

In general, the stars are believed to dominate the gravita- 
tional potential within ~Rc, with the dark matter domi 



nating outside this scale (e.g. Brighenti fc Mathews 1997 
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Figure 2. Left: Variation of the marginalized a a s a function of Rc from our data (circles), showing a strong anti-correlation. We overlay 



(stars) the lensing data of Koopmans ct al. ( |2006D , which is computed at scales smaller than ~Ro and plotted versus Re. For a strictly fair 
comparison with our work, their data may need to be shift ed s lightly to the left (given differences in the adopted photometric bands). We 
also overlay the predictions of our toy model (solid line; § [3.2| ) and an estimate of the model uncertainty (dotted lines). To indicate how 
changing the Sersic index (n) of the stellar light component affects our model, we have annotated each model line with the corresponding 
value of n. Right: Marginalized 1-cr joint confidence contours for a and M75 for each object. We overlay the predictions of our toy model, 
combined with the K-band Kormendy relation (solid line), and an estimate of the uncertainty in the model (dotted lines). 



Gerhard et al. 



2001); in p articular, we have shown this for 



the systems studied here (|H06|; |h09| [Gastaldello et al.||2007l 
Buote & Lewis "|2004|; Zappacosta et al 



Lewis et al 



2003 



2006). Therefore, the approximately powerlaw mass distri- 
butions showiriirTlg''^liidl£ate~a^^ 

tween the luminous and dark matter to produce a scale- 
free total mass distribution, at least within ~10Re. At much 
larger scales, there is evi dence that the profiles deviate from 
this s imple shape (e.g. Lewis et al 2003; Vikhlinin et al 



2006| ; [Humphrey et al.| [2006^ pastaldello et al.| |2007| ; |Ro- 

manow:;ky et al. 200£). Nevertheless, in the inner regions 



of the systems, this effect is similar to the "bulge-halo" 
conspiracy found at much smalle r scales Rc) in lensing 



and stellar dynamics studies (e.g. Treu & Koopmans 



Koopmans et al 



2004 



20061 , |00|). o appreci- 

ate that the observed powerlaw mass profiles cannot simply 
arise from the arbitrary superposition of any reasonable stel- 
lar and dark mass components. To illustrate this quantita- 
tively, we simulated a powerlaw-like mass profile and fitted 
it with a model comprising Sersi c an d dark matter compo- 
nents (as described in detail in § 3.2). For any given galaxy 
luminosity, we found that both Rc and the normalization 
of the dark matter component cannot vary by more than 
~25% from their "best" values while maintaining mean ab- 
solute fractional deviations between the m odel and data of 
less than 10% (as indicated by our fits in § 2.2). 



The origin of this conspiracy is unclear, but it is likely 
tied to the complex interaction of baryon s and dark mat- 
ter in the centres of galaxies (for example, Robertson et al 



2006 point out the importance of gas physics in maintain- 



ing the tilt of the fundamental plane, which may be related 
to this conspiracy; § i.3). Unfortunately these processes are 



very poorly understood; the predicted central dark matter 



et al 



1)86 



Gnedin et al 



densi ty cusps due to "adiabat ic contraction" (Blumenthal 
II2OO4I) n 



200q ; [Gnedin et alj p007[ ; [Putton et al.[ gOOT]), suggesting 
that other e ffects, such as dynamical friction, may be im- 
portant (e.g. [El-Zant et al.[[2C)04[; [Dutton et al.[ [2007[; [Abadi 



et al 



2009). It is unclear, therefore, that any model is so 
far robustly predicting a powerlaw total mass distribution 
arising from these effects. 

Since the data were not all fitted over the same physical 
scales, it is important to assess whether the systematic vari- 
ation in a we observe is simply an artefact of these different 
fit ranges. In Fig ^ we show the mass profiles of three rep- 
resentative systems with different Rc, arbitrarily scaled for 
clarity, and shown as a function of physical radius (left) and 
fraction of the viri al radius (right; de duced from the virial 
masses reported in Buote et al. 2007). In conjunction with 



Fig (which shows the same profiles as a function of fraction 
of Rc), it is immediately clear that throughout overlapping 
radial scales the slopes of the three mass distributions are 
all very different, so that the fit raiige is unlikely to be the 
cause of the trends observed in Fig bl 

Although all of the systems we considered are morpho- 
logically relaxed, which should eliminate objects with the 
largest departures from hydrostatic equilibrium, it is still 
important to consider whether deviations from hydrostatic 
equilibrium could be responsible for the observed trends. 
In order to affect the shape of the mass profile, any source 
of nonthermal pressure must vary radially in a finely bal- 
anced manner, given that the gas pressure profile is con- 
sistent with a powerlaw mass density profile (Fig ^ that 
varies systematical ly with the host p roper ties (FigGl). For 
the cluster A 2589, Zappacosta et al. (2006) concluded that. 



if the true mass distribution is much cuspicr than inferred 
from the X-rays (for example, an unmodified NFW plus a 
stellar mass component), no plausible source of turbulent 
or magnetic pressure could make the X-ray measurements 
imply such a fiat total density profile. 



H09 



carried out full 



ally confirmed (e.g. Humphrey et al 2006 ; ^appacosta et al" mass decompositions for three of our galaxies, finding stellar 
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mass to light ratios in excellent agreement with the predic- 
tions of stellar population synthesis models, suggesting that 
nonthermal pressure is very small (< 10-20%) within Re. 
This was further supported by the good agreement between 
the X-ray constraints on the central black hole masses, and 
those obtained from optical methods for these objects. Fur- 
thermore, n i ne of our objects were included in the sample of 
Buote et al. (2007), who showed that their dark matter halos 



inferred from X-ray mass analysis approximately exhibit the 
relationship between halo "concentration" (ci,i, ) and total 
mass predicted from numerical structure formation simula- 
tions. This supports the idea that the shapes of the mass 
profiles for these systems are not significantly in error due 
to non-hydrostatic efi'ects. 



In their survey of 15 lensing systems, Koopmans et al 



(2006) found that a ranged from ~1.8-2.4, while the Re 
spanned ~2-17 kpc. In the left hand panel of Fig ^ we 
overlay their measured a versus Re. We find that they are 
roughly consistent with the points from our X-ray sample, 
but at given Re, a is generally larger. For strict consistency 
with our adopted photometric bands, their data may need 
to be adjusted, however. Since Re in the K-band is typically 
smaller than in bluer bands (as an i llustration of this effect 



for the galaxies in the catalogue of La Barbera et al. 2008 
we find that that Re is on average ~50% larger in the R- 
band than the K -band), this would likely involve shifting the 
Koopmans et al. to the left, bringing them into better agree- 



ment with our results. By inspection, there is a slight hint in 
their data of the anti-correlation we observe between a and 
Re, although it is not statistically significant. However, since 
their fits were all restricted to regions < 0.9Rc, where the 
stellar mass dominates the potential, they are more likely to 
have been subject to local curvature in the mass distribution 
than the X-ray data were, given the better radial coverage 
of the latter. 



3.2 Decomposing the mass profiles: a toy model 

To investigate quantitatively the implications of the appar- 
ent conspiracy between dark and luminous mass components 
implied by Fig |^, we constructed a simple toy model by de- 
composing a pure pow erlaw mass distribution into an NFW 



(Navarro et al. 1997) dark matter profile plus a de Vau- 



couleurs stellar mass component (using the de projected Ser 



sic approximation of Prugniel & Simicn 1997, and fixing the 
Sersic inde x, n = 4). In our prev ious w ork ([Buote et al ' 
2007|; |H06|; |H0^; [Gastaldello et al.| ^007|; [Zappacosta et al 



2006), we have shown that decompositions of this type for 



the systems studied in the present work give parameters 
for the NFW model which are consistent with theory. The 
measured NFW concentration parameters are slightly higher 
than predicted median from N-body simulations, but this 
can be reasonably attributed to relaxed X-ray selected sys- 
tems having formed early. To perform the decomposition, 
we first constructed a set of mass "data-points" from the 
powerlaw model in 10 approximately logarithmically-spaced 
bins between ~0. 4-200 kpc, which we then fitted with the 
two-component model by minimizing the rms of the frac- 
tional residuals. The normalization and scale radius of the 
dark matter component were allowed to vary freely, as were 
Re and the normalization of the stellar model. We allowed 
the slope of the input mass distribution to vary system- 




Figure 4. Sample simulated mass profile and associated mass 
decomposition. The data-points correspond to a powerlaw mass 
distribution with o = 1.8, while the soUd fine indicates the best- 
fitting two-component (dark plus luminous matter) model. Error- 
bars, corresponding to a 10% fractional uncertainty, are shown 
to guide the eye. Also shown are the corresponding best-fitting 
dark matter (dotted line) and stellar mass (dashed line) models. 
Note how the sum of two scale-dependent models produces an 
approximately powerlaw total mass distribution over the fitted 
range. 



atically from a = 2.2 to a = 1.35. Below this limit, we 
found that the derived relationships (discussed below) be- 
came non-monotonic, which may represent the limitations of 
our simple parametrization rather than any physical insight, 
and so we prefer not to discuss this regime. While one might 
a priori expect significant degeneracies in this decomposition 
given the number of free parameters, in fact the shapes of the 
Sersic and NFW mass profiles are sufficiently different that 
we found a unique decomposition for any input powerlaw, 
at least for the range of profiles and models we investigated. 
We show an example simulated profile, and the best-fitting 
dark and stellar mass components, in Fig ^, clearly illus- 
trating how they can conspire to produce an approximately 
scale-free model. 

We overlay in Fig |2| (left panel) the best-fitting relation 
between Re and a obtained from our simple mass decom- 
positions (solid line), which can be approximated to better 
than 1% by Q = 2.31 — 0.54 log i?e. Without any fine tun- 
ing, it clearly captures the overall behaviour of the observed 
data reasonably well. Since the stellar light profiles of real 
galaxies may not be exactly de Vaucouleurs in form, we have 
also experimented with varying the Sersic index, n, from 3- 
5. Similarly, the NFW model used to parameterize the dark 
matter halo over the fitting range and may not be strictly 
correct if the baryonic and non-baryonic components inter- 
act gravitationally. There fore, we have tr i ed re placing it with 

(pool), allowing the 



the revised mass model of Navarro et al 



mass slope index to fit freely and thus provide more freedom 
in parameterizing the radial distribution of dark matter. The 
dotted lines in Fig ^ bound the range of recovered relations, 
given these changes, which clearly bracket the observed data. 
In practice, the impact of changing the stellar light profile 
dominated this estimate of our model uncertainty. 

Since our simple toy model can be freely renormalized. 
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we cannot use it alone to make any explicit predictions 
for the relationship between a and M75. One way to make 
further progress, however, is to exploit the "Ko rmendy re- 
lation " between a galaxy's luminosity and Rc (Kormendy 
1977). For a given stellar mass-to-light ratio (7*), A/75 must 



be adjusted to keep the stellar component consistent with 
the Kormendy relation while maintaining the powerlaw total 
mass distribution for a given a (and hence Re). We adopted 
a form for this relation which is appropriate for our data, 
estimated from the K-band g a laxy l umin ositie s and Rc tab- 
ulated by ICastaldello et alj ( |2007| ) and [HOel We found a 
mean relation log L/L© — 10.95 + 0.79 log i?e, with an in- 
trinsic scatter of ~0.11 dex, where L is the total luminosity 
of the galaxy in the K-band and Lq is the luminosity of the 
Sun, assuming a K-band absolute Solar magnitude of 3.41. 
We assur ned 7, = 1, app ropriate for an old stellar popula- 



tion (e.g. 



Maraston 



L , app r 

|oo§). 



In Fig Q (right panel) we overlay 



the predictions of this simple model (solid line), and an esti- 
mate of the range of uncertainty, which factors in the effects 
of varying the Sersic index of the stellar light by ±1, the 
intrinsic scatter in the Kormendy relation and a ±50% vari- 
ation in 7, (which dominates the error budget). Clearly the 
model captures the overall behaviour of the data fairly well. 

To verify that our model is self-consistent, it is impor- 
tant to determine if the parameters of the fitted dark mat- 
ter model it predicts are reasonable. Since the toy model 
cannot be used alone to set the normalization of the mass 
profile, to do this we considered the simple model, folding in 
the Kormendy relation, described above. Starting with the 
a-Mrs relation predicted by this model, we simulated cor- 
responding powerlaw mass profiles, and decomposed them 
into Sersic and dark matter components. We found mono- 
tonic relations between M^ir and a and between NFW c^ir 
and Q, allowing us to recast this as a Cvir-Mvir relation, 
which we can approximate by Cvir = C14 (M^ir/IO^^M©) 
where C14 — 10 and j3 — 0.17. This i s very close to th e 
best-fitting Cvir-Mvir relation found by |Buote et al. ( ^007| ), 
for which C14 = 9.0 ± 0.4 and /3 = 0.172 ± 0.026. Exploring 



the estimated range of uncertainty on the model linking a 
and M75, we tound C14 to lie m the range 5-15 and p from 
0.13-0.18. The total range of Alvir spanned by our predicted 
Cvir-Mvir relation (~ 12.5 < logioMvir < 15.0) is similar to 
the range spanned by the systems in our sample. It appears, 
therefore, that the parameters of the NFW component im- 
plied by our simple model are appropriate for the kinds of 
system studied in this work. 



3.3 Implications for the fundamental plane 

The existence of a conspiracy between the dark and lumi- 
nous matter to produce a powerlaw total mass profile should 
have implications for the tilt of the fundamental plane. The 
fundamental plane may be written in the form 



log Re — A log ao + Bfi + const. 



(2) 



where cro is the central velocity dispersion, jj, is the mean 
surface brightness within Rc, and A and B are proportion- 
ality constants. One way to gain insight into this relation 
is to interpret it in terms of the scalar virial theorem (e.g. 




log,„(R/kpc) 

Figure 5. The predicted ratio between the total mass and stel- 
lar mass enclosed within Rc from our toy model (7m), shown 
as a function of Re (solid line). The dot-dash lines indicate 
our estimate of the uncertainty on the model, and the dashed 
line indicates the best-fitting linear approximation in the range 
logi?c =0.2—1.0. The d ata-points are the observ ed K-band M/L 
ratios from the data of La Barbera et al. 

I daoosl ) (see text). Note 
the good agreement between the data and the toy model. 



Binney & Tremaine 2008). We can write this in the form: 
log Re = 2 log ao - log (J^^ - log 7 + 5.47 



Here Re is in kpc, ao is in km s~^, i'i/2 is (approximately) 
the light within Rc (in Solar units) and 7 = 7,7m is the mass 
to light (M/L) ratio in Solar units, where 7* is the stellar 
M/L ratio and 7^ the ratio of total mass to stellar mass 
within Rc. The value of the additive constant on the right 
hand side is useful only in setting the absolute normalization 
of 7, and does not affect our conclusions. For convenien ce we 
chose t o det ermine it from the theoretical relation of Woli 
( 2005: ) ^, siii ce this implies a mean 7* of ~1 for the 



et al 



La Barbera et al. 



(2008) galaxy sample (see below), which 



is consistent with the expectations for an old stellar popu- 
lation. Alternativ e , emp irical relations have been found by 
Cappcllari et al.| (^00e| , assuming mass follows light) and 
' p007| , , 



Bolton et al. ( 2007| , assuming a singular isothermal sphere 
mass distribution) , which would imp ly lower 7* , by as much 
as ~40% (for the Cappellari et al.| relation). Nevertheless, 



for our work the absolute value of 7* is unimportant, and so 
this choice does not affect our conclusions. 

In the K-band, assuming the absolute Solar magnitude 
is 3.41 and neglecting cosmological effects, we can write 
log(Li /2/'7ri?e) ~ 16.00 — 0.4/i, and hence 



Jl/2 

log Re ^2 log (To + OAfi - log 7 - 10.53 



(3) 



If structural non-homology does not significantly affect this 
equation, the tilt of the fundamental plane arises through 
the dependence of 7 on the other properties of the galaxy. 
Eliminating ctq between Eqns |^ and ^ and rearranging, we 



^ We replaced the total velocity dispersion in their relation with 
the central velocity dispersion, which is only strictly valid for a 



flat velocity dispersion profile. 
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obtain: 
log 7 = 



('213 \ 
0.4j /i + const. 

(0.33 ± 0.04) log + (0.007 ± 0.011)/i + const. 
(0.33 ± 0.04) log i?e + const. (4) 



where we have used the "log q- q fit" values for A and B from 
the recent K-band analysis of La Barbera et al, (2008). The 
^ term is consistent with zero, within the error, suggesting 
that the tilt of the fundamental plane can arise if 7 depends 
only on Rc, which is exactly the prediction of the toy model, 
as we discuss below. Nevertheless, since the ^ term may be 
non-negligible, to obtain a more accurate assessment of how 
7 varies with Re, we computed log 7 directly from the data 
of La Barbera et al. by using Eqn ^ We grouped these data 

into a gpripg nf narrrnnf ^'''^'j^, ^"'^^^''i ParVi nf ^nrViirVi wp 



calcula Jed the mean -y and its error (by applying the central Barbera et al 



that the observed behaviour could be explained in terms 
of the predicted curvature. Indeed, there is a hint that the 
observed 7 versus Rc relation begins to deviate from our 
simple predictions at Rc<; 20 kpc (Fig which may arise 
due to deviations from structural homology among BCGs, 
or it may represent the limitations of our simple model. 

To make a more quantitative comparison between 
the traditional fundamental plane and our toy model, we 
adopted a linear approximation for the toy model relation- 
ship between 7„i and Rc (log7m — —0.11 + 0.371og7?e; this 
is accurate to about ~0.03 dex over the range 0.2 < log Re < 
1.0, as shown in Fig ^, assuming 7, — 1). Substituting this 
into Eqn |^ (assuming 7* is constant) and rearranging into 
the form of of Eqn|, we found A=1.46 and B=0.29. These 
are in good agreement with recent K-band measurements 
of the fundamental plane (A~l 4-1. 5 and B~0. 30-0.33: H 



limit theorem to the scatter of the data), as shown in Fig ^ 
We can use our toy model not only to predict the rela- 
tionship between Re and a, but also the dark matter fraction 
within Re (and hence 7m), as is immediately clear from in- 
spection of Fig ^. Since the dark matter fraction does not 
depend on the overall normalization of the mass profile, we 
do not need to fold in any other relations (such as the Kor- 
mendy relation). Furthermore, since we find that Re varies 
monotonically with a (Fig it is possible to recast this 
as a relation between 7™ and Rc. We show the resulting 
relation (roughly log 7™ — —0.03 log i?e + 0.32 (log_Re)^) in 
Fig 1^ along with an estimate of the uncertainty that arises 
from from varying th e Sersic index by ±1, or adopting the 



Pahre et al 



2008 



Jun & Im 



2008 



Mobasher et al. 



199c 



1995). Repeating the analysis, but allowing the 



Sersic index to vary by ±1, o r usin g the alternative dark 



matter model of Navarro et al. (2004), we obtained A in the 
range ~1.2-1.6 and B~0.23-0.31. Thus, assuming only a 
powerlaw total mass distribution, we can clearly reproduce 
the shape of the K-band fundamental plane, without any 
fine-tuning of the model parameters. We therefore speculate 
that the establishment of a powerlaw total mass distribu- 
tion is a fundamental feature of galaxy formation and the 
principal factor in determining the tilt of the fundamental 
plane. Since the relative importance of gas-rich and "dry" 
merging for the assem bly of present-day giant ellipticals is 



still under debate ( e.g. Hopkins et al 



Navarro et al 



(pOM) model for the dark matter profile. It [Nipoti et alj |2009|), the requirement of producing a nearly 



2008 



Naab et al. 



200c 



is immediately cle ar that there is go o d agr eement between 



the model and the La Barbera et al. (2008) data, implying 
the mean 7* ~ 1, as noted above. Since the stellar popula- 
tions of early-type galaxies (and hence 7«) do not strongl y 



vary with Re (at least for Rc?J 1 kpc: Graves et al. 2009) 
this suggests that assuming the ubiquity of approximately 
powerlaw density profiles (supported by Fig |l|) is sufficient 
(albeit not necessary) to explain the tilt of the fundamen- 
tal plane. While it is possible that not all giant early-type 
galaxies have density profiles which so closely resemble a 
powerlaw (although we are not aware of any convincing such 
counter-examples found by X-ray methods), to lie on the 
fundamental plane. Fig |^ indicates that their 7™ must still 
vary with Re in a very similar way to the predictions of the 
toy model. Nevertheless, since it is actually the flattening 
of the density proflle with increasing Re that is responsi- 
ble for the 7m-Ro relation predicted by the toy model, it is 
plausible that similar flattening of the mass distributions in 
a system without a powerlaw-like mass distribution might 
produce approximately the same trend. 

One intriguing prediction of our toy model is that the 
relationship between log7m and log7?e has some curvature, 
which would also imply curvature in the fundamental plane. 
There is, in fact, some evidence that brightest cluster galax- 
ies (BCGs) may not lie on the same fundamental plane as 



powerlaw total mass distribution in the core of low-redshift 
systems will therefore be a key test of these models. 
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